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A new model of localized highly frustrated ferromagnetism is presented: kagome spin ice. By use of 
analytical and Monte Carlo calculations its massive groundstate entropy is evaluated. Monte Carlo 
calculations are also used to explore the phases in the presence of further-neighbor interactions and 
applied fields. The importance of thermal spin fluctuations and their ability to stabilize a variety 
of magnetic structures is clearly manifested. Remarkably, some of these phases are only partially 
ordered and present both disordered and ordered sublattices. 



PACS number: 75.10.-b, 75.40.-s ,75.25.+z 
I. INTRODUCTION 



Geometric frustration refers to the inability of a mag- 
netic system to minimize its individual exchange interac- 
tions, as a simple consequence of the exchange geometry. 
For many years it was thought to be limited to antifer- 
romagnets, with much research focusing on the kagome 
and pyrochlore antiferromagnets, 1-5 until a phase termed 
'spin ice' was discovered by Harris and coworkers. 6-8 In 
their system, which is based on the pyrochlore lattice, the 
large number of degenerate ground states that character- 
ize frustrated magnetism is a consequence of particular 
Ising anisotropics and ferromagnetic nearest-neighbor ex- 
change. Mapping the spin directions onto the disordered 
hydrogen bonds of Pauling's cubic ice model, their system 
was shown to remain disordered even at T=0. However, 
if other terms are present, such as dipolar interactions 
a transition to long-range order has been found due to 
collective zero-energy flips of loops of spins. 9 Despite the 
unique properties of pyrochlore spin ice, nothing is known 
about 2-dimensional (2-d) spin ice models or more gener- 
ally about the properties of frustrated 2-d ferromagnetic 
systems. In contrast, the ice models were widely inves- 
tigated in two dimensions, especially on a square lattice 
a version of which was solved exactly by Lieb. 10,11 Ex- 
tensions and generalizations of the model and to other 
lattices soon followed. 11,12 

In this article we examine the single layer kagome fer- 
romagnet with in-plane Ising anisotropies and identify a 
new spin ice system — kagome spin ice. Using analytic 
and classical Monte Carlo calculations we demonstrate 
the massive groundstate degeneracy of this ice-like phase, 
which greatly exceeds that of the pyrochlore spin ice. 
Furthermore, we show how its degeneracy is raised by 
additional parameters in the Hamiltonian, namely sec- 
ond neighbor interactions or an applied magnetic field, 
and that these lead to the stabilization of a number of 
exotic magnetic phases. 




FIG. 1. a) The 6 degenerate spin configurations for a single 
triangle of the kagome spin ice model with the Ising axes 
pointing to the center of each triangle. The Bethe lattice b) 
and the kagome lattice c) of vertex sharing triangles. The 
nearest-neighbor and second neighbor interactions, Ji and J2 
respectively, of the kagome lattice are defined separately in 



II. THE ICE RULES 

These local conditions define whether a spin configu- 
ration lies within the ground state manifold. For the py- 
rochlore system, 6-8 where the Ising axes point towards 
the centre of the basic tetrahedral motif, the ice rules are 
that 2 spins enter and 2 spins leave each tetrahedron. In 
the case of the 2-d kagome lattice, where the Ising axes 
now point towards the centre of the triangular motif, the 
ice rules are that either 2 spins enter and one leaves each 
triangle, or the converse. Inspection quickly reveals that 
there are 6 possible ground state configurations associ- 
ated with each triangular motif. These are demonstrated 
in Fig. la). Extension of the motifs in 2 dimensions shows 
that kagome spin ice possesses an extensive degeneracy, 
the residual entropy of which is calculated in Section III. 
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III. RESIDUAL ENTROPY A LA PAULING 



IV. MONTE CARLO CALCULATIONS 



The residual entropy of cubic ice, and thus so of py- 
rochlorc spin ice, has been evaluated by Pauling 13 as 
S p = While this value is very close to the exact 

one evaluated numerically, 14 we note that in fact Paul- 
ing's estimation is that of a Bethe lattice of tetrahedra: 
starting from an arbitrarily chosen 'central' tetrahedron, 
6 of the 2 4 = 16 possible spin states correspond to spin 
ice configurations. As the 4 neighboring tetrahedra each 
have 1 spin fixed due to sharing with the central tetra- 
hedron, they each have only 3 possible ice states, from 
a complete manifold of 2 3 = 8 states. The same argu- 
ment can then be used for all tetrahedra as one moves out 
through the Bethe lattice. Thus, among the 2 N possible 
states (N being the total number of sites), only a frac- 
tion (|) T (Nt being the number of tetrahedra, i.e. 
= N/2) satisfy the ice rule on all tetrahedra, which 

leads to an entropy: S p = kln(2 N x (f)^) = ^Zn3/2. 
In the pyrochlore lattice, closed loops have to be taken 
into account: for one tetrahedron with 2 spins fixed by 
the neighboring tetrahedra, the number of acceptable ice 
states is either 1 or 2 (depending on the configuration of 
the 2 fixed spins), among 4 possible configurations for the 
2 remaining spins; this leads again to an average reduc- 
tion factor of | , while if 3 spins are fixed this factor is 
larger: \. Thus for a lattice of tetrahedra, the number of 
ice states should be slightly larger than the Pauling value 
Sp , as has been found numerically S n — LOIS',,. 14 . As 
previously demonstrated for more general ice models 12 , 
Pauling's estimate is a lower bound on the exact value of 
the entropy. 

The estimation of the residual entropy for a Bethe lat- 
tice of triangles (Fig. lb) can be made following the same 
principles: for each triangle with 1 fixed spin, the number 
of ice states is equal to 3, while the number of possible 
states is 2 2 = 4. Thus, if Nt is now the number of trian- 
gles, the total number of ice states is 2 N x (|) T , (with 
Nt = ^), leading to an entropy S b = ^Zn9/2. If we 
again consider the effect of the closed loops present in the 
kagome lattice (Fig. lc), we find that if 2 spins are fixed 
for a given triangle, the number of acceptable states for 
the third spin is either 1 or 2, among 2 possible states. 
Assuming that the states with 2 spins fixed are equiprob- 
able, an average factor | is again obtained. Although this 
procedure ignores any correlations beyond the first neigh- 
bors, this simplified argument suggests that closed loops 
should not be expected to change that much the entropy 
in the kagome lattice and the exact value should be very 
close to Sb- Our numerical estimations show clearly that 
this is the case. In fact this spin ice kagome system can 
be mapped exactly on the Ising kagome antiferromagnet, 
in a similar way as the pyrochlore system. 15,16 Conse- 
quently, the value that we found for the entropy is the 
same that was found for disordered ground state of the 
AF Ising kagome system. 17,18 



Our classical simulations were made on a single 
kagome plane of 50 x 50 unit cells, each containing 3 
spins of modulus 5 = 1, i.e. 7 500 spins in total. The 
nearest- and second neighbor exchange integrals, J\ and 
J2, were used to couple the spins according to the path- 
ways shown in Fig. Id. A Kirkpatrik 20 cooling scheme 
was used with each step in temperature involving a reduc- 
tion from Ti to T2 according to the equation T2 = 0.9T\. 
At each temperature 5 000 cycles site -1 were used for 
equilibration and 5 000 cycles site -1 for the calculation 
of the thermodynamic properties. The starting temper- 
ature for all the simulations was T = 10 1 J\\S 2 . The spin 
maps presented are representative of the average over the 
Monte Carlo cycles performed at each temperature that 
are within the specified range. 

The Hamiltonian can be defined as: 

H = -J2 ■ $3 + h a- (9Vb) ( l ) 

i,j<i i 

where the anisotropy term Ha reads ^2 Ki [Si ■ Xi) . 

i 

Ki < is the single-ion anisotropy constant that forces 
the local spin Si to align along the local axis Xi and H 
is the applied field. In order to create Ising spins with 
the different local orientations shown in Fig. 1, the trans- 
verse spin fluctuations were completely frozen out, which 
corresponds to setting = —00. 
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FIG. 2. (a) The specific heat and the (b) magnetic entropy 
of kagome spin ice (•) and the related antiferromagnetic sys- 
tem (o) as a function of temperature. Due to their differing 
scales, the comparison between the specific heat of the two 
systems is shown as an insert. The limiting values of their 
magnetic entropies, Smag, as T — » 00 are shown. 

V. THERMODYNAMIC PROPERTIES AND SPIN 
CORRELATIONS OF KAGOME SPIN ICE 

The specific heat of kagome spin ice shows a broad 
maximum at T/J\S 2 = 1. Above this temperature, in 
the range 1 < T/J1S 2 < 2, the slow approach of the 
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magnetic entropy to its saturation level reveals the pres- 
ence of significant short-ranged order correlations. The 
residual entropy of the spin ice state is clearly discernible 
in the magnetic entropy (Fig. 2) as the shortfall in the 
high temperature limit with respect to the ln2 expected 
for Ising spins as T — > oo. If the limiting value of the in- 
tegrated entropy is taken as 0.1948 (35), we find that the 
residual entropy of the spin ice phase is S mag {T = 0) = 
S mag (T = oo)-/ °° §dT = Zn2-0.1948 = 0.4982 (35), a 
value that is in good agreement with the residual entropy 
of Sb = |/n9/2 = 0.5014 calculated for the Bethe lattice. 

As a contrast, the entropy of the antiferromagnctic sys- 
tem ( Ji < 0) with the same Ising anisotropics as the spin 
ice phase is also shown in Fig. 2. While a maximum in 
the specific heat is also seen at Tc — \J\S 2 \ (Fig. 2a), 
the magnetic entropy differs sharply from the kagome 
spin ice in that it involves a recovery of the complete 
S = Inl value expected for the formation of long-range 
order. The long-range order that results is the so-called 
q=0 spin structure familiar in kagome antiferromagnets 2 
and is drawn in Fig. 3. 




FIG. 3. The 'q=0' ground state of the kagome antiferro- 
magnet. 

Equal time correlations < Si ■ Sj > between distant 
spins Sj and Sj, averaged over the Monte Carlo cycles 
at equilibrium, were computed at different temperatures. 
Fig. 4 show the thermal variations obtained for first and 
second spin neighbors. An increase in absolute values is 
observed as the temperature is decreased, down to T ~ 
0.5 1 Ji|S 2 at which the correlations saturate. 
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FIG. 4. A semi-log plot of the first (•) and second neighbor 
(□) spin-spin correlations in kagome spin ice as a function of 
temperature. 

The correlations at zero temperature in the spin ice 
phase are evaluated analytically for the Bethe lattice of 
triangles for all neighbors as < Si ■ Sj >b— (-jp) cos &ij > 
where m is the number of bonds that form the path join- 
ing the spins Si and Sj, and #y is the angle difference 
9j — Oi between the orientations of the local anisotropy 
axes Xj and Xj for these spins. Numerical values for dif- 
ferent neighbors Sp, S 7 , . . . of a given spin S a as spec- 
ified in Fig. 5 are reported in Table I, and compared 
to the correlations for corresponding neighbors in the 
kagome lattice, that were computed at the lowest tem- 
perature (T ~ 0.0004 JiS 2 ) of the Monte Carlo simula- 
tions. Agreement is remarkable, except for the S a — S$ 
spin pairs. This is to be expected as in the kagome lattice 
these neighbors belong to the same hexagon and are thus 
connected by two paths with exactly the same number of 
bonds, which leads to correlations that are larger by a 
factor of ~ 2. In a similar way the other distant spins 
in the kagome lattice, such as the second neighbors, are 
also not connected by only one path, however as these 
other paths have a larger number of bonds and the con- 
comitant corrections decay as a power law in the number 
of bonds, they give rise to only minor adjustments. 




FIG. 5. Neighborhood of a given spin S a in the Bethe and 
the kagome lattice of triangles. 
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TABLE I. Spin-spin correlations in the Bethe and kagome 
lattice of triangles between a given spin S a and its neigh- 
bors as specified in Fig. 5 (n th : order of neighborhood in the 
kagome lattice). 



As J2/J1 is increased, the intermediate phase becomes 
stable over a wider temperature range; when J2/J1 = 
—0.4 this phase occurs at temperatures from O^JiS 2 to 
0.6 JiS 2 . 



a) b) 




FIG. 7. Averaged Monte Carlo spin maps of the structures 
for J 2 < when a) \J 2 \S 2 <T < 1.5|J 2 |5 2 b) T « \J 2 \S 2 . 



Open circles represent disordered spins. 



VI. ANTIFERROMAGNETIC SECOND 
NEIGHBOR INTERACTIONS 

These were examined with the values of J2/J1 = -0.4, - 
0.2 and -0.1. In all cases 3 features are seen in the specific 
heat (Fig. 6). The first at T ~ J1S 2 involves the forma- 
tion of short-ranged spin ice-type correlations. Below this 
temperature, at T ~ 1.5|Ja|5 2 , the additional antiferro- 
magnetic exchange stabilizes a partially ordered struc- 
ture that is made up of a triangular motif of clockwise 
and counter-clockwise hexagons of spins (Fig. 7a). The 
mean magnitude of the moments on the sublattice that 
separate these hexagons is zero. The reason for some of 
the sites failing to order is made clear from examination 
of the lowest maximum at T ~ | J2IS' 2 : this involves the 
formation of an ordered magnetic structure which has a 
unit cell that is 3 x larger along each of the axes than the 
paramagnetic cell (Fig. 7b). The complexity of this struc- 
ture is evidenced from the 2 types of hexagon: the clock- 
wise hexagons (and the time-reversed counter-clockwise) 
are separated by those that have 3 pairs of parallel spins. 
As the latter involve adjacent spins with antiferromag- 
netic 120° orientations this ground state configuration 
is frustrated. Re-examination of the intermediate phase 
(|J 2 |S 2 < T < l.5\J 2 \S 2 ) shows that the fluctuations 
have selectively annulled the mean moment on sites that 
would otherwise frustrate the lattice. 
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FIG. 6. Specific heat for the Ising kagome ferromagnet with 
second neighbor interactions: J2 < 0. 



VII. FERROMAGNETIC SECOND NEIGHBOR 
INTERACTIONS 

As the lattice of second neighbors is also a kagome lat- 
tice, ferromagnetic J2 still result in a highly frustrated 
spin system and if the moments related by J\ and J2 are 
of uniform magnitude 1/3 of the moments lie on sites 
of zero mean-field (Fig. 8). It is not surprising, there- 
fore, that extensive fluctuations are still present at the 
lowest temperatures for all the values of J 2 > studied 
(J2/J1 = 0.4, 0.2, and 0.1), in sharp contrast with the 
ground state when J 2 < 0. If these second neighbor inter- 
actions present a minor contribution to the Hamiltonian, 
e.g. J2/J1 = 0.1, short-ranged ice-like correlations still 
form and these give rise to a broad feature in the specific 
heat at T ~ J\S 2 (Fig. 9). However, in contrast with 
the J2 < phases, only one maximum is seen at lower 
temperature, T ~ J2S 2 . Remarkably, this is not associ- 
ated with a phase transition and no long-range order is 
seen down to the lowest temperature: only a tendency 
is observed that involves the condensation of a partially 
ordered phase with the propagation vector k=(0 0) and a 
non-zero ferromagnetic component (Fig. 8). These corre- 
lations are strongest when J2 is weakest. As J2 increases, 
this tendency dwindles out and the two maxima in the 
specific heat merge. What remains are strongly fluctuat- 
ing phases with correlations that are less and less remi- 
niscent of the partially ordered phase. 21 
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FIG. 8. While no long-range order occurs, tendencies to- 
wards order with particular structures are seen. These spin 
maps are presented for J 2 > when a) T < J 2 S 2 b) 
T << J2S 2 . Open circles represent disordered spins. 
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FIG. 9. Specific heat for the Ising kagome ferromagnet with 
second neighbor interactions: J2 > 0. 
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FIG. 11. The reduced magnetization per site as a function 
of applied field along the a direction (a) in units of the molec- 
ular field JiS and (b) after rescaling the applied field with 
temperature, i.e. plotting M/Mo vs. [iBBg /kT , a universal 
curve is obtained. The temperatures of the field sweeps were 
(o) 0.254, (□) 0.051, (•) 0.010, and (■) 0.002 JiS 2 . The lines 
provides guides to the eye. 



VIII. EFFECTS OF MAGNETIC FIELD 

The effects of an externally applied magnetic field on 
the spin ice phase were studied with fields of up to 
0.269 J±S applied along the characteristic directions of 
the kagome lattice, i.e. a and a* (Fig. 10). These calcu- 
lations were made at a temperature of 0.025 J\S 2 after 
the system had beencooled in zero applied field. 




FIG 



a) H //a b) H //a* 

10. The high field spin configurations when the ap- 



plied field is along the a and a* directions. Open circles rep- 
resent disordered spins, i.e. those that have equal probability 
of being in either direction. 



With the field parallel to the a direction a continu- 
ous increase in an ordered component is observed. No- 
tably, it involves only 2/3 of the sites (Fig. 10a). This 
demonstrates another point of similarity between the 2-d 
kagome spin ice system and that of the closely related 
3-d pyrochlore lattice: the spin components induced by 
a field along some of the crystallographic directions are 
incapable of choosing a unique ground state and signif- 
icant degeneracy is retained. A common result is the 
formation of only partially-ordered spin configurations. 

With the field parallel to the a* direction a continuous 
increase of order with the k=(0 0) structure (Fig. 10b) is 
again found. Trivially we find that the ordering is pref- 
erential for the sites with their local anisotropy parallel 
to the field direction. The ordered component of all the 
sites saturates in a field of 0.054 JiS. 
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FIG. 12. The reduced magnetization per site as a func- 
tion of applied field along the a* direction (a) in units of the 
molecular field JiS and (b) after rescaling the applied field 
with temperature, i.e. plotting M/Mo vs. fisBg/kT, a uni- 
versal curve is obtained. The temperatures of the field sweeps 
were (o) 0.254, (□) 0.051, (•) 0.010, and (■) 0.002 J1S 2 . The 
lines provides guides to the eye. 
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Inspection of the high- field ordered structures suggests 
that even an infinitesimal field applied along either of the 
characteristic directions would induce an ordered compo- 
nent. It is therefore a surprise that the initial dependence 
of the reduced magnetization on the field is not sharper. 
Our simulations show that this weaker dependence is a 
result of thermal fluctuations, rather than being a conse- 
quence of the finite lattice size. This is demonstrated by 
the collapse of the reduced magnetization as a function of 
the Zeeman energy divided by the thermal energy, onto 
the same curve at all temperatures studied (see Fig. 11 
and Fig. 12). 

IX. DISCUSSION 

The thermodynamic properties of kagome spin ice, and 
those of phases in the presence of second neighbor in- 
teractions and applied fields, manifest clearly the rich- 
ness of frustrated physics once thought to be the do- 
main only of antiferromagncts. In particular, the im- 
portance of spin fluctuations is central and their pres- 
ence leads to the stabilization of a variety of magnetic 
structures with both disordered and ordered sublattices. 
While such partially ordered phases have already found 
to be stable in frustrated lattices when the local moments 
are close to an instability 22 or in systems with contrived 
further-neighbor exchange, 23 here they arc stable with 
well defined magnetic moments and are a simple conse- 
quence of the strong frustration associated with uniax- 
ial anisotropy. In real systems 24 additional terms such 
as quantum fluctuations and dipolar interactions will be 
present in the Hamiltonian. While further work is re- 
quired to determine how these will alter the spin ice phase 
and the partially ordered phases presented here, these ef- 
fects have already been shown to be strongly influential 
on their low-temperature physics. A review of how Ising 
models are influenced by the introduction of quantum dy- 
namics that arise from the inclusion of an XY exchange 
component is presented in Ref. 26. An interesting limit 
for their inclusion is the XY antiferromagnetic on a py- 
rochlore lattice model presented by Champion et al. 27 , in 
which spin-wave analysis indicated that zero-point quan- 
tum fluctuations stabilize the formation of an ordered 
moment from the a highly degenerate continuous mani- 
fold. 

The observation that the ordered Neel state proposed 28 
for the vanadium jarosite (H30)V3(S04)2(OH)6 is a 
highly symmetric member of the spin-ice manifold sug- 
gests that spin ice-like correlations will be present above, 
but close to, its ordering temperature. The system is 
clearly not a true realization of kagome spin ice as the 
presence of a transition to Neel order requires that either 
additional terms in the exchange Hamiltonian that we 
have not considered here have raised the degeneracy of 
the ice-like manifold, or that the single-ion anisotropies it 
possesses are of a magnitude that quantum fluctuations 



are important. Unfortunately, the single- ion terms in 
the jarosite are difficult to estimate because the relative 
contributions of axial (-DS" 2 ) and rhombic (i?(5 2 - S* 2 )) 
terms determined from magnetic susceptibility data 29 
taken from KGa2.96Vo.o4(S04)2(OH) 6 are ambiguous, as 
other pairs (D, E) parameters may fit them equally well. 
It is still, however, reasonable that the short-range cor- 
relations present above the ordering temperature will be 
characteristic of the most important terms in the Hamil- 
tonian, and thus of the spin ice phase. 

In conclusion, we have presented a new model of local- 
ized highly frustrated ferromagnetism: kagome spin ice. 
If residual entropy is used as a measure of frustration, our 
analytical and numerical calculations demonstrate that 
kagome spin ice is more highly frustrated than the 3-d 
counterpart, pyrochlore spin ice. It is therefore the most 
highly geometrically frustrated "ferromagnetic" ground 
state yet studied. Furthermore, second neighbor inter- 
actions within the kagome plane are shown to have the 
salient effect of stabilizing structures that feature both 
ordered and disordered sublattices. 
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gold Laboratories, University College London, 20 Gordon 
Street, London, WC1H OAJ, UK. 
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